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Abstract 

We present new developments of the simple model of the central engine of GRB, proposed recently. 
The model is based on minimal assumptions: some rotating compact relativistic object at the center 
and stable perturbations of its rotating gravitational field, described by Teukolsky Master Equation. 
We show that using nonstandard polynomial solutions to the angular Teukolsky equation we can 
describe the formation of coUimated jets of various forms. Appearance of imaginary part of the 
superradiance-like frequency is established for the first time for pure vacuum black hole jet solutions 
of Teukolsky equation. 

1 Introduction 

Gamma-ray bursts (GRB) are among the most powerful astrophysical objects in the Universe. They are 
considered as highly coUimated {9 jet ^ 2° — 5°) explosions on cosmic distances emitting huge amounts 
of energy in very short periods of time seconds) [1, 2, 3, 4]. Despite the increasing amount of 
observational data, the theory of GRB is still far from being clear. One of the major problems is the 
lack of understanding of the nature of the so called central engine - its physical nature and the process 
involved in the emission of huge amounts of energy (~ 10^^ — 10^^ erg). 

Unexpected feature discovered by the mission SWIFT is the existence of flares - a surprising evidence 
of late time activity of the central engine [3, 4, •'3]. There are clear indications that in some cases the flares 
are produced by energy injection by the central engine. Although number of theories and models have 
been suggested to explain the light curves of GRB, the process however remains a mystery that cannot 
be solved without a good model of the central engine. 

The current focus of GRB physics is on the propagation of the emitted matter that can explain the 
observed light curves. The most used is the fireball model (for details see [1]). That model, however, 
cannot explain the observed flares. Another basic model is the cannonball model of long GRBs [6] . 

The central engine - the physical object producing the GRB - before SWIFT epoch usually was 
considered to be a Kerr [7] black hole (BH) - a hypothetical result of the death of a massive star for the 
long GRB (Tgo > 2s) (see [■]], [n]) or as part of a binary merger of compact objects (BH-BH, NS-BH, 
or NS-NS; where NS=neutron star) for short GRB (Tgo < 2s). The first hypotheses seems to be hardly 
compatible with the observed flares, produced via energy injection by the central engine. The second 
hypotheses was recently refuted by the existing detectors of gravitational waves [9] . The latest analysis 
shows that the short and long GRB may have a similar central engine (except for its duration) [" ]. 

It was believed that one of the possible models for such engine could be a Kerr black hole (KBH) in 
super-radiant mode [11, 12] - the wave analogue of the Penrose process. The late-time evolution of a 
perturbation of a Kerr metrics is governed by the quasi-normal modes (QNM) - oscillation with complex 
frequencies that are determined only by the BH parameters (mass, charge and angular momentum) 
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[13, 11]. Those frequencies are obtained solving the Teukolsky Master Equations with certain boundary 
conditions. According to the standard theory, supperradiance occurs in the process of scattering of 
exterior waves on KBH for rea/ frequencies lo < iOcriUcai = am/2Mr^ (m > - integer, > being the 
event horizon radii) [111]. 

Another process which was currently related with the GRB central engine is the Blandford-Znajek 
process [14] based on electromagnetic extraction of energy from KBH. Sometimes it is considered as an 
electromagnetic analog to Penrous process. The realization of this idea has its own achievements and 
problems (see in [14]). 

However, the Penrose process seems to be not efficient enough to explain the huge amount of promptly 
radiated energy in relativistic jets (~ 10**^ — 10**^ erg) [15], as shown long time ago by Wald (see in [12]). 

In a series of talks [IG], we have presented a mathematical description of relativistic jets that can 
help to penetrate into the mysterious physics behind the central engine. A simple model, based on novel 
solutions of the angular Teukolsky equation, however, not regular, but singular ones [16, 17, 19] was 
proposed. They seem to describe in natural way the most important feature of the relativistic jets - their 
coUimation. In this paper, we would like to announce the recent results of our numerical simulations. 
The main new result is the appearance of a complex critical frequency, which play the role of the well 
known supperradiance ones. Its imaginary part is of the same order of amplitude as the real part and 
yields an exponential abatement of the superradiance-like emission in jets, created by KBH. 

2 A Toy Model of Central Engine 

Our simplified model was already presented in [Ifi, 17]. Thus, we will omit the details and give only a 
summary of the basic ideas behind it. In brief, the jets observed in GRBs require a rotating object - KBH 
or a compact massive matter object, to produce them. In the jet's problem we can use the Kerr metric to 
describe the gravitational field of that object at least in a very good approximation [17], since the visible 
jets are formed at distances about 20-100 event-horizon-radii [15]. At such distances one practically is 
not able to distinguish the exterior field of KBH from the exterior field of rotating objects of complete 
different kinds (see for example [ : ] and the references therein) . 

As a result, the only working way to get information about the real nature of the central engine is 
to study jet's spectra. The different objects yield different frequencies of perturbations of space-time 
geometry, because of the different boundary conditions on their surface (See the articles by Fiziev in [13] 
and [17]). Thus, measuring the real jet's frequencies one can get indisputable evidences about the actual 
nature of the central engine. 

In the present article we are probing only the KBH model for generation of jets. The results for 
other possible models will be discussed elsewhere. Due to the well known "no hair" theorems, the non- 
charged KBH solutions of Einstein equations, which potentially may be of astrophysical interest, depend 
on extremely small number of parameters - the mass and the angular momentum. The unique boundary 
conditions on the event horizon yield unique robust spectra, defined only by the mass and the angular 
momentum of the KBH. These spectra differ essentially from the spectra of relativistic jets, generated 
by central engine of any other nature. Thus the KBH-based model of relativistic jets seems to be easily 
recognizable from observational point of view, at least up to the disguise effects by the KBH environment. 

Using the Teukolsky Master Equation (TME) we acquire the linearized perturbation of the Kerr 
metric (see [11]). From there, we follow the procedure established by Teukolsky for separation of the 
the variables in the TME using the substitution: <i> = e('^*+™'*)'S'(0)i?(r) where m = 0,±1,... and 
Lu — LOR + iuji is a complex frequency (with ujj > 0). 

The exact solution of radial equation for the radial function R(r): 

A^R{r) + 2 {s + l){r- M)^R{r)+ ( ~ '^''^ ~ ^ - 4isuj r - x] R{r) = (1) 
av' ar \ A J 

can be expressed in terms of confluent Heun functions [16, 17, 18, 19]). 



2 



The most important for GRB physics seem to be electromagnetic waves with s = — 1. In this case 
two independent exact solutions in outer domain are: 



-+1 



Ri{r) = Cie""^'' {r~r^) [r - r^) HeunC [ a, P,-f,d,T],-- ^) (2) 

and 



i?2 (r) = Cge"''^'^ (r- r_^)^' (r- r_)' HeunC (^a, -/3, 7, 5, ry, -^J— ^ j 



(3) 



/ \ „ 2i(w(a^+r^)+am) 2 i(w (a^ + r_ ^)+am) 

where: a = (^r^ — r_j tj, p = ^ _^ 1,7= F~^r ' 

(5 = — 2i(r^ — r_)aj(l — i{r_ + r^) oj), 

?7 ^ — 4a;2r^4 ^_ ^4^^ _ gw^r ) r_^3 + (l - Aaujm - 2iJ-a^ - 2A) [r^"^ + r_2) + 

(^r^-r j L 

(4 io) r_ - 8iuj + 4:A- iuj^a"^ ~ 2) r_ - 4a^ (m + w a)^ 

For the angular function S{6), as an exact solution the angular equation: 



o . (fn + su) 
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(auju) + 2aujsu + s + 



5/™ = (4) 



we use polynomial solutions in terms of Heun polynomials ([16, 17, 18, 19]). In equation (4) we are using 
the variable u = cos 9. Putting f2 = aw we have explicit expressions: 

= e±" (cos ((?/2))l^-"l (sin (0/2))-l^+'"l x (5) 



HeunC ( ±4ri, | s-to|, | s+m|, -4f2s, h 2 s - fi^ - A - s, cos 



2 ^ 
2 

and 

41)^ „^(0) = e±" ^''^^^ (cos (0/2))l^-"l (sin (0/2))-l^+'"l x (6) 
HeunC ( T4ri, | s+m\, \s-m\,-Afl s, 2^s ~ fl'^ ~ A - s,sin^ - 

The polynomial solutions of TME describe one way waves and are most suitable for modeling of the 
relativistic jets [17]. Posing the polynomial condition on the angular solutions, we obtain for the one way 
electromagnetic waves on Kerr background the explicit formula 

^=-i,m(w) = - 2 m±2 Vfi^ + fi TO. (7) 

This simple form of the separation constant A for polynomial solutions is the most significant mathemat- 
ical advantage of our model of jets, generated by electromagnetic perturbations of Kerr metric. 

The KBH boundary condition for the radial equation ca be obtained using the following assumptions: 

1. On the horizon we allow only incoming in the horizon waves. Then we obtain the different solution 
working in each interval of frequencies: a) for m — 0, only R2; b) for m > if luu G (— cJcnO) and 
R2 on the outside; for m < i?i when luu G (OjCJcr) and R2 outside. 

2. On infinity we allow only outgoing waves. In general, the function i? is a linear combination of an 
ingoing (i?<_) and an outgoing (R^) wave: R = R^ + R^ with some constants C<_ and C_^. In 
order to have only outgoing waves, we need to have — 0. This equation defines the spectral condition 
for the frequency lu. The main mathematical problem is that the explicit form of the constant is 
not known. Therefore to solve the spectral problem we use the around way, proposed by Fiziev in [ ]: 
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The straightforward check shows that for solutions (2), (3) we have Hm = in the special direction 
arg{r) — 3t:/2 — arg{uj) in complex plane C,.- Taking lim in this direction, we obviously obtain the 
spectral condition for uj in the form 

= hm = 0. (8) 

3 Numerical Results 

The numerical evaluations with confluent Heun functions are generally too complicated for a number 
of reasons. At present the only software package able to deal with them is Maple. Unfortunately, the 
existing versions of Maple package require too much time for numerical evaluation, even on modern fast 
computers with large amount of memory. Additional problem is that the procedure calculating those 
functions is not working well in the whole complex plane and special attention must be paid to the 
regions where the procedure becomes unstable. 



3.1 Visual Resemblance of the Solutions and Real Jets 

We already reported that plotting the solutions of the angular equations with relation (7), we obtain 
figures that resemble jets (see [IG, 18]). Illustrations of the jet features of the model can be seen from 
the animations of such oscillating solutions on our site [ ]. 

An interesting comparison with Nature can be obtained looking at the picture of the discovered by 
NASA's Spitzer Space Telescope "tornado-like" object Herbig-Haro 49/50, created from the shock waves 
of powerful protostellar jet hitting the circum-stellar medium. "More observations should help us to 
unravel its mysterious nature" [21]. Without any doubts in this case the jet is not related with BH. This 
one, as well as many other real observations prove the presence of BH to be not necessary for generation 
of jets. The collimated jets of GRB are another example of a possible application of our simple model. 
There the main problem still remains the physical nature of the central engine, too. Clearly, if possible, 
a common model of jets of different scales is most desirable. 




Figure 1: Some of the shapes of jets that our model can represent. Compare it with the picture of 
Herbig-Haro 49/50 (see [:'1]) 



3.2 Complex Super-Radiance-like Jet Modes 

The study of the solutions of radial Teukolsky equation (1) is much more complicated from computational 
point of view. Even after fixing the relation (7) between A and lo, solving the KBH boundary condition (8) 
for radial equation is not a straight-forward process, because of the analytical and numerical problems 
connected with the evaluation of the confluent Heun function. Our method consisted in plotting and 
examining the spectral condition in the complex plane C^^ to find points that resemble roots of the 
transcendental equation (8). Then we tested those points with more precise root-finding algorithms. 
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Thus wc found for the case s — —1, 2M — 1, a/M — 0.99, \m\ — 1 two roots: 

cji,2 = ±0.8676087 + i 0.1236275. 

For the case m = we obtained another pair of roots 

tJi,2 = ±0.006394 + i 0.1354325801. 

To confirm those roots, we plotted the studied function in small regions around the zeroes and obtained 
the expected conus-like form of the modulus of the function in the complex domain and the complete 
rotation of it's phases in angle 2tt around it's simple zeroes. 

Interesting in these roots is that the absolute value of their real part is precisely the known critical 
frequency tUcr = am/2AIr^,m > that is connected with the superradiance phenomena (for a recent 
review see the last reference in [1 '■]). To the best of our knowledge this is the first time that an imaginary 
part of this quantity is discovered solving BH boundary conditions in pure vacuum (i.e. without any 
mirrors, etc). In the case of jet's solution we have a critical phenomenon around the complex value of the 
frequency. It's explicit manifestation is the change of the role of solutions (2) and (3) outside the even 
horizon of KBH when \ujr\ crosses the critical value am/2Mr^. The complexity of the critical frequency 
shows an essential difference between our jets-from-KBH solutions and standard QNM, obtained using 
regular solutions of the angular equation (4) [ ]. It is also important that the real and the imaginary 
parts of the critical frequency are with the same magnitude. 

The lack of zeros with < in our numerical studies can be both due to the stability of the model, 
or to numerical problems in that area. The study of the zeroes is a field of intensive research by the team 
and more conclusive results will be presented somewhere else. 

4 Conclusion 

Our simple toy model of central engine seems to be able to produce some of the basic features observed 
in GRB. The Teukolsky Master Equation with the correct boundary conditions is fundamental enough 
to account for all types of GRB. The polynomial angular solutions of TME show a jet-like structure. In 
the radial equation, the essential assumptions we used is that the imaginary part of the frequency with 
BH boundary conditions (only entering waves on the horizon and only going to infinity waves) should 
be positive. This provides stability of the solutions in direction of time-future infinity and indicates an 
explosion in direction of time-past infinity. The roots we found numerically agree with our assumption. 
Interesting new result is the complex critical frequency w^'^* — ujn+iuj (loj ^ ujji) of " superradiance" that 
appeared, showing that jet's to-be-superradiance modes decrease exponentially in time in the direction 
of the future and blow up in direction of the past. 
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